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ABSTRACT 

Solar flares produce hard X-ray emission of which the photon spectrum is of- 
ten represented by a combination of thermal and power-law distributions. How- 
ever, the estimates of the number and total energy of non-thermal electrons are 
sensitive to the determination of the power-law cutoff energy. Here we revisit 
an 'above-the-loop' coronal source observed by RHESSI on 2007 December 31 
and show that a kappa distribution model can also be used to fit its spectrum. 
Because the kappa distribution has a Maxwellian-like core in addition to the 
high-energy power-law tail, the emission measure and temperature of the instan- 
taneous electrons can be derived without assuming the cutoff energy. Moreover, 
the non-thermal fractions of electron number/energy densities can be uniquely 
estimated because they are functions of the power-law index only. With the 
kappa distribution model, we estimated that the total electron density of the 
coronal source region was ~ 2.4 x 10 10 cm" 3 . We also estimated without as- 
suming the source volume that a moderate fraction (~ 20%) of electrons in the 
source region was non-thermal and carried ~ 52% of the total electron energy. 
The temperature was 28 MK, and the power-law index 6 of the electron density 
distribution was -4.3. These results are compared to the conventional power-law 
models with and without a thermal core component. 

Subject headings: Sun: flares — Sun: particle emission — Sun: X-rays, gamma 
rays 
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Introduction 



A solar flare is an explosive energy release phenomenon on the Sun and accelerates a large 



number of electrons up to tens of MeV (e.g. lBrownlll97ll ; iLin fc Hudson! 1 19 76c iMiller et al. 



19971 ; iHolman et al.l 120031 ) . To diagnose accelerated electrons, the hard X-ray (HXR) obser- 
vations of electron bremsstrahlung emission have been used. 

In general, the spatially integrated HXR photon spectrum exhibits a relatively flat , non- 
thermal tail in addition to an intense and steep thermal component (ILin et al.lll98lf). Al- 
thoug h a model with multiple temperatures can often fit the entire spectra (jEmslie fc Brown 
19801 ). the non-thermal tail can typically be described as a power- law or double power-law 
that connects with the thermal component typically in the 15-30 keV range. 

When viewed as an image, the intense thermal emission is dominated by an arcade 
loop structure whereas the less intense but the high energy tail of the HXR emission is usu- 
ally d etected at chromospheric footpoint of the loop (e.g. 



Hoyng et al. 



1983). The non-therm al tail can also originate from the corona (e.g. Frost fc Dennis! 



Palmer &: Smerd 



1994 



1981 



Brown et al. 



19721). and the source is somet imes located 'above-the-loop' (e.g. lMasuda et al 



1971 



Krucker et all 120081 : llshikawa et all 120111 ). 



A caveat of studying the non-thermal HXR emission is that, the thermal emission from 
a loop is so bright that it masks spectral features of the non-thermal sources (in either 
footpoint or 'above-the-loop') especially in the lower energy range. Therefore, it is difficult 
to clarify how far the power-law spectrum extends in the lower energy direction. As such, a 
low energy cutoff E c of the power-law has been considered, typically in the 15-30 keV range, 
to estimate the number and total energy of non-thermal electrons in the source, although 
the estimates can be sensitive to the choice of E c . 

Thus, efforts have been made to under stand properties of HXR emission around E c 



(e.g. IHolman fc Benkal Il992l ; ISui et al.l 120051 ). In particular, since the launch of RHESSI, 
it has been successfully shown that a range of values for E c fit the data equally well and 
that the highest value of E c that still fits the dat a can be used to derive the lower-limit 
for the non-thermal numb e r and energy densit ies ( IHolman et al.l l2003t lEmslie et al.l 12004 ; 
Saint-Hilaire &: Benzl 120051 ; iKontar et al.l 120081 ) . As for the physical meaning of the low- 
energy cutoff, it was argued that the cutoff represents the critical velocity above which elec- 
trons run away and are freely accelerated by the reconnection electric field ( IHolman &: Benka 
19921 ). If a sharp cutoff existed, however, plasma instabilities wo uld lead to flattening of the 
distribution around and below the cutoff energy (as reviewed by lHolman et al.ll201ll ). 



In fact, a theoretical study pointed out that the non-thermal electron distribution could 
seamlessly merge into a thermal distribution, removing the need for a low-energy cutoff 
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( lEmslid 120031 ) . Moreover, in-situ observations of electrons in the planetary and interplan- 
etary space often show that the higher energy tail of a thermal core component smoothly 
exte nds into a power-law distribution. Examples can be found at the Earth's bow shock 



(e.g. 



Gosling et al.lll989l ) and the magnetotail reconnection (e.g. l0ieroset et al.ll2002l ). 



In some cases of in-situ observations, the kappa distribution model (jVasyliunasI Il968[ ) 

has been used to represent the entire elect ron distribution because it is a composite of 

Maxwellian-like core and a power-law t ail (e.g. Christon et al.|[l988 . 1989 . 1991 ; Onsager fc Thomsen 



199ll : lLeubner!l2004j ; 



Imada et al. 



20 



l|) . While the kappa distribution was first introduced as 



an empirical model ( Vasyliunas] 1968 ). recent theoretical and computational studies have sug- 
gested that self-consistent formation of electron kappa distribution is po ssible by the beam- 
plasma interactions which invo l ve the Langmuir /ion- sound turbulence ( jYoon. Rhee &: Ryu 
20061 ; iRhee. Ryu fc Yoonll2006l ; iRyu et al.l 120071 ). The origin of the kappa dist ribution has 
also been discussed in terms of Gibbsian theory (jTreumann fc Jaroschekl 120081 ) and Tsallis 
Statistical Mechanics ( Livadiotis fc McComas 201 ll. and references therein). From the solar 
physics point of view, IKasparova fc Karlickyl (120091 ) has already suggested that the kappa 
distribution may also be useful for interpretations of solar HXR sources. They reported 
that a kappa distribution fits the spectrum of a coronal loop-top source but fits less well the 
spatially integrated spectrum (coronal and footpoint sources). 



The purpose of this paper is to complement the work of IKasparova fc Karlickyl (120091 ) 
by examining the kappa distributi on model in a recently reported RHESSI event of 2007 
December 31 (IKrucker et al.ll2010l ). We studied this event because an unusually intense 
HXR emission was detected from 'above-the-loop' coronal source. We show that the spatially 
integrated HXR spectrum can be fitted by not only a combination of thermal and thin-target 
power-law distributions but also a combination of thermal and kappa distributions. The 
introduction of the core distribution in the non-thermal source via the kappa distribution 
enables us to estimate the number and energy densities without assuming the cutoff energy. 



2. Kappa Distribution 

The isotropic, three-dimensional (3D) form of the kappa distribution function f K {v) (s 



cm 6 ) is written as 

Jk[V) (tt/^ 2 ) 3 / 2 r(« - 1 /2) V + k9 2 J 1 ' 

where v is the particle speed, k is the power-law index, 9 is the most probable particle 
speed, T is the Gamma function and N K is the number density. The coefficient is deter- 
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mined so that J f(v)d 3 v = N K . If k is sufficiently large, the distribution approaches a single 
Maxwellian distribution. The most probable energy is E mp = (\/2)m8 2 at which the dif- 
ferential flux (=(v 2 /m)f(v)) becomes maximum. However, the temperature is defined as 
k B T K = {l/2)m6 2 [k/(k — 3/2)] so that the average energy of particles can be expressed as 
-Eavg = (3/2)fcsT K . Note that E avg = (3/2)A;bTm for the isotropic 3D Maxwellian distribution 
Im{v) with UbTm = (l/2)m,Vt h where v t h is the thermal speed. 

By using the kappa temperature k B T K and introducing particle energy E = (l/2)mv 2 , 
we can convert f K (v) into the density distribution F K (E) (cm -3 keV -1 ) as 

2N K y/E T(k + 1) 



V«(k B T K )z (k _ 3 /2) 3/2 r(R _ 1/2) 



E 



- -(k+1) 

X 1+ k B T K (K-?>/2)\ (2) 
so that J F K (E)dE = N K . The thi n-target formula of this expre ssion has been incorporated 



into the Solar Soft Ware (SSW) by lKasparova fc Karlickyl (120091 ) and can be used as OSPEX 
fitting function f _thin_kappa.pro. 

An example of F K (E) is plotted in Figure H^a). The Maxwellian distribution Fm{E) 
with kTu = E mp = kT K [(n — 3/2)/ k) is also plotted for comparison. Here, the density Nm 
of the Maxwellian distribution F M (E) has been adjusted so that F K (E mp ) = F M (E mp ). Such 
N M is derived as 



n k r(«-i/2) V k. 

Then, the difference between F K (E) and the adjusted Fm{E) represents the non-thermal 
particles, as indicated by the shaded region. A slight difference remains in the lower energy 
range (E < E mp ) but the total difference in this range is negligible compared to the total 
difference in the higher energy range (E > E mp ). Therefore, the ratio Rn of the non- 
thermal electron density to the total electron density in a source can be approximated by 
Rn = 1 — Nm/N^- The ratio R £ of the non-thermal electron energy to the total electron 
energy can also be calculated using -E avg - Both R^ and R £ are functions of k only and are 
shown in Figure [0(b). 

For the example case of k = 4 (i.e. Figure QJa)), the non-thermal electrons constitute 
~20% of the total electrons and such non-thermal electrons carry ~50% of the total electron 
energy. Note that we do not need to assume cutoff energy nor source volume to estimate 
the values. Also, these ratios are much larger than what is generally assumed in an electron 
beam model where beam density is much less than the ambient density ('diluted beam'). 
Presence of a kappa distribution in a coronal source region may indicate that a significant 
number of electrons are locally accelerated. 
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3. Analysis 

To test the kappa distribution model, we performed imaging spectroscopy for a partially 
disk-occulted s olar f lare of 2007 December 31 observed by RHESSI. Following the work by 



Krucker et all fcOlOf ). our focus goes to the time of HXR peak flux 00:47:42-00:47:50 UT. 



Figure |2] shows the spatial structure of the HXR sources during the peak time in 8 
different energy ranges. There were mainly two separate sources. The northern source at 
(X,Y)~(-980,-150) arcsec was dominated by a low energy (<10 keV) X-ray emission (Figure 
|2{a,b)) whereas the southern source at (X,Y)~(-970,-165) arcsec was dominated by a high 
energy (>20 keV) X-ray emission (Figure |2fg,h)). The northern source in the 6-8 and 8-10 
keV ranges was so bright that the flux of the southern source should be less than what can 
be calculated in, for example, the blue polygons of Figure [2]Ja,b). Conversely, the fluxes in 
the red polygons of Figure E^g,h) would be the upper limits of the northern sources in the 
20-25 and 25-30 keV ranges. 

Note that, soon after the HXR peak time, the main thermal loop of the southern source 
was identified along the limb on the western side of the southern source. Thus, the northern 
source has been considered as a separate thermal source, although the precise relation be- 



tween the northern and southern sources remains unclear (IKrucker et al.l 120101 ) . Below, we 
focus on the spectral features of the sources rather than a possible relationship between the 
two sources. 

Figure [2]^ c-f) show the details of the HXR sources in the intermediate energy ranges. In 
the 10-12 and 12-14 keV ranges (Figure [2b and[2H, respectively), both northern and southern 
sources appear together, indicating that both sources had comparable fluxes. In the 14-17 
and 17-20 keV ranges (Figure [2fe and [2f, respectively), only one source can be identified 
somewhere between the northern and southern sources and it is not clear which of the two 
sources it belongs to. Because of such unclear nature of the source structure, we took the sum 
of the fluxes in the dashed polygons and consider it as the upper limit of both northern and 
southern sources. The sum, however, is essentially the same as the values in the spatially 
integrated spectrum shown below and so we do not use the data from Figures EJje) and 
|2]T) in the following analysis. Then, assuming that the red and blue polygons represent the 
northern and southern sources, respectively, we calculated the total flux within each polygon 
to be compared with spectrum models. 

It is to be noted that, while we chose the polygons so that the double sources in Figure 
|2^c) can be separated, another choice of boundary indicated by the white line in Figure [2^d) 
could also be used to better separate the double sources in Figure EJ^d). We found that such 
modification to the polygons leads to less than 30% flux change in all energy ranges. Thus, 
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we will use this number as the uncertainty of the measured fluxes. 

Figure [3] shows the result of the imaging spectroscopy. The light-red and light-blue 
squares indicate the fluxes from red and blue polygons in Figure |2J and the fluxes are 
compared to the spatially integrated photon spectrum (histograms) as well as four different 
model distributions (colored curves). The steeper and flatter components are evident in the 
integrated spectrum and we will use thermal, power-law and kappa distributions to represent 
these two components. To represent the non-thermal tail, Models A and B use the power-law 
distributions whereas Models C and C use the kappa distribution. 

For the power-law distribution fits used in Models A and B, we used a formula that calcu- 
lates thin-target Bremsstrahlung X-ray spectrum from a power-law electron flux density dis- 
tribution (cm~ 2 s~ 1 keV~ 1 ). This formula is implemented as f _thin.pro in the SSW/OSPEX 
software and contains three free parameters a 55 , 5 FD and E c . a 55 (10 55 electrons cm _2 s _1 ) 
is the normalization factor and is a product of the number density of plasma ions, the flux 
density of non-thermal electrons, and the volume of the radiating source region. 5fd is the 
power-law index of the flux density distribution. Throughout this paper, however, we use 
the power-law index 5 of the number density distribution when comparing different models. 
The two indices can be converted to one another by 5 = <5fd + 0.5. E c (keV) is the low 
energy cutoff of the power-law distribution. 

For the kappa distribution fits used in Models C and C, we used the SSW/OSPEX 
procedure f _thin_kappa . pro but it does not contain line emissions. We added line emissions 
in the analysis and imposed that the same values of emission measure (EM) and temperature 
(T) are used in the kappa distribution and line emissions. It is to be noted that the assumed 
values for the fits were derived for a Maxwellian distribution, and the temperature inferred 
under the assumption that electron distribution is Maxwellian may be an overestimate o f 



the actual temperature of a distribution with non-thermal tail ( lOwocki fc Scudderl Il983l ). 
As for the spectral index, the program f _thin_kappa . pro assumes a kappa distribution for 
the electron density and the spectral index k can be converted to 5 by 5 = k + 0.5. 

Figure |3] Model A uses the thermal (red curve) and power-law (blue curve) distributions 
for the steeper and flatter components, respectively. The light-red and light-blue squares 
are consistent with the red and blue curves, respectively, indicating that the northern source 
was producing the steeper component whereas the southern source was producing the flatter 
component. Note that a range of values for the low-energy cutoff E c fit the data equally well. 
We found that the highest E c that still fits the data is 16 keV with \ 2 — 1-0 ( n °t shown). 
However, this model underestimates the flux in the 10-12 keV range by 46^5 % °f what 
was measured in the image (light-blue square). Then, we imposed E c to be in the range 
11 < E c < 13 keV in order to look for a solution consistent with the imaging spectroscopy 
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result. We found that, as shown in Figure [3] Model A, E c = 12 keV would best fit the data, 
although x 2 became relatively large (x 2 = 1-3). 

Figure [3] Model B uses the thermal (red curve) and combination of thermal and power- 
law (blue curve) distributions for the steeper and flatter components, respectively. This 
is basically the same as Model A, but an additional thermal distribution is introduced to 
account fo r the possible c ore co mponent of the southern source (blue curve). A similar model 
is used by ICaspi fc Linl (120101 ) . This set of distributions can also fit the data nicely (x 2 = 
0.9) including the 10-12 keV range, suggesting that a core distribution could have existed in 
the southern source. Again, a range of values for the low-energy cutoff E c fit the data equally 
well and we chose the highest E c that still fits the data (E c = 35 keV). Note that the model 
curves (generated to fit the integrated photon spectrum) have a higher energy resolution so 
that the modeled line emissions partially exceed the fluxes obtained from images in the 6-8 
and 8-10 keV ranges. We confirmed that the modeled values averaged over the same energy 
ranges are consistent with the values from the imaging spectroscopy. In Model B, the thermal 
core component has a relatively large temperature of ~52 MK a nd this is comparabl e to the 
temperature ~44 MK of a super- hot coronal source reported by lCaspi fc Linl (120101 ). 

Figure [3] Model C uses the thermal (red curve) and kappa distributions (blue curve) for 
the steeper and flatter components, respectively. Although the reduced \ 2 is slightly larger 
(x 2 = 1.2), this set of models can also represent the data nicely including the 10-12 keV 
range. Again, the modeled line emissions partially exceed the fluxes obtained from images 
in the 6-8 and 8-10 keV ranges, but we confirmed that the modeled values averaged over 
the same energy ranges are consistent with the imaging spectroscopy. The number of free 
parameters, 5, is still the same as that of Model A (thermal + power-law) whereas Model B 
(thermal + thermal + power-law) needed 7 parameters to have a thermal core distribution 
in the flatter component (blue curve). 

Figure [3] Model C uses the kappa distributions for both steeper and flatter components. 
It is evident that the non-thermal tail of the steeper component is still below the upper limits 
(light-red arrows). The slope is quite soft (k ~12), however, indicating that the steeper 
component (northern source) was mostly thermal. Note also that the core temperature of 
the steeper component 10 MK is reduced by 50% compared to the temperature from the 
thermal fit, 21 MK. While such low temperature of the steeper component is still comparable 
to the temperature 15 MK measured by GOES during the same interval (00:47:42-00:47:50), 
the emission measure obtained by the fit was unrealistically high, 2.5xl0 49 cm -3 , and it is 
an order of magnitude higher than what was measured by GOES, 1.5 x 10 48 cm -3 . Therefore, 
this set of two kappa distributions is not favorable for representing the observation. 



To better visualize the above comparisons, Figure 0] uses the model curves in the 5-31 
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keV range and takes the ratio of natter component (blue curves in Figure |3]) to steeper 
component (red curves in Figure [3]). The ratios are compared to the flux ratios from images 
(black squares with error bars). The modeled values are averaged over the energy ranges of 
the images. It is evident that all models are consistent with the imaging spectroscopy. The 
upper limits in the < 10 keV energy range indicate that, in this energy range, X-ray emission 
from thermal plasma in the southern source could not have been detected by RHESSI if its 
flux is less than 65% of the flux from the northern source. 

Figure O examines the time variations of the fitted parameters for the flatter component 
(i.e. southern source). While the peak time spectra in Figures [3] were taken during the 8 
s interval of 00:47:42-00:47:50 UT, we fitted the data every 4 s in Figure [51 The spectral 



indic es S from all models show soft-hard-soft variation as already reported by (IKrucker et al. 



20101 ). The low-energy cutoff E c of Model B (power-law with thermal core) is systematically 
larger than that of Model A (power-law without thermal core) because of the presence of the 
thermal core distribution in the flatter component (southern source). As for the parameters 
of the core distribution, the emission measure EM and the temperature T of Model B 
(thermal-l- power-law) are systematically lower and higher, respectively, than those of Model 
C (kappa). Therefore, Model B suggests presence of a superhot thermal core distribution in 
the flatter component (southern source) whereas Model C suggests presence of larger number 
of non-thermal electrons. The reduced \ 2 fluctuated around ~ 1, indicating that Models A, 
B and C fitted the data fairly well. The averages in the shown interval (00:47:15 - 00:49:10 
UT) are Model A: 1.22, Model B: 0.88 and Model C: 0.94. 



Discussion 



Let us now discuss implications of the results based on Models A, B and C for the peak 
flux interval (00:47:42 - 00:47:50 UT). We will particularly discuss non-thermal fractions 
of electron number/energy densities in the southern source. The estimated non-thermal 
fractions are summarized in Table 3. 

Model A (thermal + power-law) assumes that the southern source ('above-the-loop' 
coronal source) contains a negligible amount of thermal electrons and uses the power-law 
with no thermal core to represent the flatter component. To estimate the number density of 
non- thermal electrons ('i nstantaneous' density), we assumed a source volume V = 8 x 10 26 
cm 3 ( IKrucker et al.ll2010l ) and applied the formula by lLinl (119741 ) to the power-law part of the 
photon spectrum. To estimate the number density of thermal electrons ('ambient density'), 
we assumed that the ambient envi ronment should be si milar to that of the nearby thermal 
source. Following the derivation by IKrucker et al.l ( 120101 ). the ratio N nt /N t ^ can be expressed 
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as 

where A^ ppcr = 8xl0 9 cm" 3 is the upper limit of the ambient density and we used 7 = 4.4, 
flux at 50 keV of 0.16 ph s _1 cm -2 keV -1 and the low-energy cutoff E c = 12 keV. If we 
use N t h = jV" t T per , we obtain Rn = N nt /N tot ~ 0.05. If we use the best estimate of N t h = 



2xl0 9 cm 3 (jKrucker et al.ll2010i ). we obtain ~ 0.44. Similarly, the ratio e nt /e t h can be 
expressed as 

£^0.47^rV(^r S (5) 



wher e we assumed that the temperature of the ambient plasma is 22 MK ( jKrucker et al. 



20101 ) . Then, if we use N t h = , we obtain R £ ~ 0.32. If we use the best estimate 

of Nth = 2xl0 9 cm -3 , we obtain R £ ~ 0.88. Because of this large fraction of non-thermal 
electrons, Model A implies that the non-thermal electrons are not simply a tail on the 
thermal distribution and that electrons are accelerated locally in the southern source (i.e. 
'above-the-loop' coronal source). However, Model A resulted in a relatively large \ 2 (~l-3 
at the flux peak time and ~1.2 on average). Therefore, we explored other possible models 
as described below. 

Model B (thermal + thermal + power-law) assumes that the southern source ('above- 
the-loop' coronal source) contains a significant amount of thermal electrons and uses the 
power-law with a hot (52 MK) thermal core to represent the flatter component. To estimate 
the non-thermal fraction of electron number density, we can use the obtained normalization 
factor a 55 (~ 0.3 x 10 55 cm _1 s _1 ) because it is actually a product of the number density 
of plasma ions, the flux density of non-thermal electrons, and the volume of the radiating 
source region. Using t he plasma density 4 .5 x 10 9 cm -3 (Table [2]) and an assumed source 



volume of 8xl0 26 cm 3 (jKrucker et al.ll2010i ). the electron flux density is estimated to be 8 x 
10 17 cm~ 2 s _1 . Based on the low-energy cutoff energy E c = 35 keV, the mean speed of the 
accelerated electrons can be estimated on the order of 10 10 cm/s and the number density of 
the non-thermal electrons in the southern X-ray source is estimated to be 8 x 10 7 cm -3 . This 
is only 2% of the thermal electron density. 

Such small fraction of non-thermal electron density is consistent with an electron beam 
scenario in which electrons are accelerated above the hot flare loops and stream through the 
source region to produce a super-hot thermal plasma. Model B is also consistent with our 
imaging spectroscopy result especially in the 10-12 keV range, justifying our assumption that 
a thermal core distribution may have existed in the southern X-ray source ('above-the-loop' 
coronal source). In fact, Model B gives the least x 2 (= 0.9) and, as such, seems to be a 
plausible model. 
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A caveat is that, because it contains an additional thermal distribution, the number 
of free parameters, 7, is relatively larger compared to 5 in Models A and C. In general, a 
larger number of free parameters contributes to decreasing the reduced \ 2 - Therefore, Model 
B may have resulted in the lower x 2 partly because of the smaller number of parameters, 
although the number of free parameters alone does not explain the \ 2 difference. It is also 
to be noted that a range of values for E c fit the data equally well and the highest value of 
E c that still fits the data has been used in Model B. As such, the above estimation of the 
fraction of non-thermal electrons, 2%, should only be considered as a lower-limit. 

In both power-law Models A and B, we needed to assume the source volume V as well 
as the low-energy cutoff E c to estimate the fraction of non-thermal electrons. Note that 
we can only obtain a very rough estimate of V and the estimated fraction of non-thermal 
electrons can be sensitive to the choice of E c . Then, we consider the kappa distribution as an 
alternative because it contains a thermal core component that seamlessly extend to a power- 
law distribution. The kappa distribution allows us to estimate the fraction of non-thermal 
electron density /energy without assuming the source volume (Section [2]). 



In Model C (thermal + kappa), the ratio of non-thermal electron density to the total 

o.oi 

0.01' 



electron density in the southern source was = 0.20iooi> an d the ratio of non-thermal 



energy to the total electron energy in the southern source was R £ = 0.521q'q2- ft is to be 
emphasized that the non-thermal fractions of number/energy densities have been derived less 
ambiguously than Model A and that the result does not invoke the possibility of non-thermal 
electrons outnumbering thermal electrons. On the other hand, the derived estimate of 
= 0.20lo oi i s m uch larger than ~ 0.02 of Model B as derived above. This implies that 
not all electrons are thermalized in the southern ('above-the-loop') source region and that 
there may have been local acceleration of electrons in this region. 



If we assume the source volume of 8xl0 26 cm 3 ( jKrucker et al.l 120101 ). the total density 
of the southern source N to t can be estimated as (1.0 ± 0.9) xlO 10 cm -3 . Within the error 
range, the estimated density is consistent with the upper limit of 8 x 10 9 cm -3 derived by 



Krucker et al.l (l2010f ). The estimated density is also consistent with what was estimated in 
Model B, (4.5 ± 2.2) xlO 9 cm" 3 (Table W). 

As for temperature, the 'above-the-loop' region (i.e. southern source) had a temper- 
ature of 28 ± 9 MK. This is ~ 1.3 times larger than the temperature 21 ± 0.4 MK of 
the nearby thermal sour ce but is ~ . 6 tim es the temperature of the super-hot component 



discussed in Model B or lCaspi fc Linl (120101 ) . We speculate that the released magnetic field 
energy was converted to both thermal and non-thermal energies of electrons but a significant 
fraction (R £ ~ 0.5) went to non-thermal electrons so that the temperature did not increase 
considerably. 
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As for the effective plasma beta (3, iKrucker et al.l (120 10[ ) estimated the magnetic field 



strength \B\ to be 30 - 50 G and derived (3 between ~0.005 and ~0.02 for a pre-flare plasma 
with the density 2 x 10 9 cm ~ 3 and the temperature 2 MK. They argued that the preflare 
thermal plasma could be replaced with non-thermal electrons (the low-energy cutoff at 16 
keV) so that the effective plasma beta becomes /3 ~ 1. If we use the density ~10 10 cm -3 
and the kappa temperature ~ 28 MK as derived from Model C, the effective plasma beta 
falls between ~ 1 and ~ 3. 

From the spectral fit, we obtained k ~3.8 which leads to 5 ~4.3 for the density dis- 
tribution F(E) oc E~ s (see Equation (2)). This is somewhat softer compared to 5 ~3.9 
derived from Model A (thermal + power-law) and 5 ~3.8 derived from Model B (thermal + 
thermal + power-law). The k- value is not too large, however, and the kappa distribution is 
far from a single Maxwellian. Thus, an electron acceleration theory still needs to reproduce 
a power-law tail for this event. Note again that the power- law index S can be converted from 
<5fd of the flux density distribution used in f_thin.pro (Models A & B; 5 = <5fd + 0.5) and 
from k of the number density distribution used in f _thin_kappa . pro (Models C & D; 5 = k, 
+ 0.5). 

It is to be mentioned that Model C also has caveats and disadvantages. First, x 2 was 
relatively larger (x 2 ~ 1.2 at the peak flux interval 00:47:42 - 00:47:50. The relatively large 
X 2 may imply that the kappa distribution is still not the best functional form to represent the 
HXR spectrum from the southern source. However, the reduced x 2 averaged over the larger 
interval 00:47:15 - 00:49:10 was smaller ~ 0.9 and this is comparable to the average x 2 ~ 
0.9 of Model B. The time variation of x 2 was also similar between Models B and C (Figure 
[5]). Thus, Model C is as acceptable as Model B in terms of spectral fitting. Second, the 
spectral index at the time of peak flux S ~4.3 is even m ore inconsistent with that deduced 
from the radio observations (5 ~3.4; Krucker et ali J2OI0I )) than that deduced from the other 



models (5 ~3.9 in Model A and 6 ~3.8 in Model B). However, the radio emission represents 
electrons with energies larger than ~ 100 keV whereas our analysis was made in the < 100 
keV range. To further understand if these disadvantages are common in other solar flare 
events, it is important to test the kappa distribution in a larger number of solar flare events. 



5. Conclusion 

The kappa distribution does not require a low-energy cutoff E c to represent non-thermal 
electrons, and the thermal core component can seamlessly extend to a power-law distri- 
bution. Furthermore, the non-thermal fractions of electron number/energy densities can 
be uniquely estimated because they are functions of the power-law index k only. While 



-12- 



Kasparova & Karlickyl (120091 ) applied the kappa distribution to loop-top coronal sources, 



we examined the kappa distribution model in an unusually bright 'above-the-loop' coronal 



source obtained on 2007 December 31 (IKrucker et al.ll2010l ). For comparison, we also exam- 



ined the conventional power-law models with and without a thermal core distribution in the 
source. 

Model A, the power-law with no thermal core component, was consistent with the imag- 
ing spectroscopy result when we chose E c = 12 keV, although the reduced \ 2 was relatively 
large (~1.2 on average). This model implies that non-thermal electrons can outnumber 
thermal electrons. 

Model B, the power-law combined with a super-hot (52 MK) thermal core component, 
could fit the observed spectrum well (x 2 ~ 0.9 on average) and was consistent with the 
imaging spectroscopy. This model implies that at least 2% of the source electrons carried 
non-thermal energies. 

However, both Models A and B require a low-energy cutoff E c to represent the non- 
thermal tail, and the estimates of the electron number/energy densities can be sensitive to 
the choice of E c . Furthermore, a source volume V had to be assumed for the estimates but 
we can only obtain a rough estimate of V. 

Thus, we examined Model C (the kappa distribution model). We found that it can 
fit the observed spectrum well (x 2 ~ 0.9 on average) and is consistent with the imaging 
spectroscopy result. Without assuming the source volume V and the lower-energy cutoff E c , 
we estimated that a moderate fraction (20%) of the source electrons had non-thermal energies 
and carried 52% of the total electron energy in the 'above-the-loop' coronal source region. 
The temperature was 28 MK and the power-law index of the electron density distribution 
was -4.3. It would be important to examine a larger number of events in order to verify the 
generality of the kappa distribution model. 
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Fig. 1. — Properties of the kappa distribution F K (E) compared with the adjusted Maxwellian 
distribution Fm(E): (a) An example spectrum for N K = 1.0, k = 4.0, T K —E mp [K/(K — 3/2)], 
N M = 0.8 and T M =E mp where E mp = 3.0 keV = 35 MK. (b) The density fraction i? N and 
the energy fraction R £ of non-thermal electrons, plotted as functions of k. 
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Fig. 2.— Hard X-ray sources during 00:47:42-00:47:50 UT of 2007 December 31. The images 
are constructed for eight different energy ranges by the CLEAN algorithm. Sub collimators 
3-8 are used with natural weighting. The color-code shows the photon flux normalized by 
the peak flux within each image. The peak fluxes were (a) 22, (b) 2.1, (c) 0.52, (d) 0.17, 
(e) 0.12, (f) 0.05, (g) 0.046 and (h) 0.016 [photons cm" 2 s" 1 arcsec" 2 ]. The blue and red 
polygons indicate the regions in which the total fluxes were calculated for the spectra in 
Figure [U The dashed polygons in panels (e) and (f) emphasize that the data from these two 
energy ranges are not used in the analysis. We confirmed that the conclusion of this paper 
is not sensitive to the choice of subcollimators, construction algorithm and the size/shape of 
the polygons. The black solid curve gives the location of the photospheric limb. 
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Fig. 3. — Comparison of the observed spatially integrated photon spectrum (histograms) 
with modeled distributions (solid curves) as well as the imaging spectroscopy result (light-red 
and light-blue squares for the northern and southern sources, respectively). Four different 
sets of models were used to fit the observed spectrum: (A) thermal + power-law, (B) thermal 
+ thermal + power-law, (C) thermal + kappa and (C) kappa + kappa. See Tables [1] and |2] for 
the resultant parameter values of the steeper (red) and flatter (blue) component, respectively. 
The gray curve is the background. The histograms in the lower panels are the residuals of 
the fit in units of the standard deviation derived from photon statistics. The bremsstrahlung 
emission formula of power-law, i.e. f_thin.pro is used for Models A and B. This program 
is implemented as double power- law but we fixed the break energy at 1000 keV to make it 
a single power-law. All kappa distributions take into account the line emissions. We used 
f_vth.pro with the lines-only option and combined it with f _thin_kappa . pro through a 
wrapper routine. 
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Table 1. Fit result of the steeper component (colored red in Figure [3]) 



fit parameters 


Model A 


Model B 


Model C 


Model C 




thermal 


thermal 


thermal 


kappa 


emission measure EM, cm" 3 


(1.2 ± 0.1) X 10 48 


(8.6 ± 1.0) x 10 47 


(4.8 ±0.2) x 10 47 


(2.5 ± 1.0) x 10 49 


total density 51 Ntot, cm -3 


(3.9 ± 1.3) x 10 10 


(3.3 ± 1.1) x 10 10 


(2.4 ±0.6) X 10 10 


(1.8 ± 1.1) x 10 11 


temperature b T, MK 


17 ±0.5 


18 ±0.6 


21 ±0.4 


10 ±0.7 


power-law index c S 








12 ± 0.8 



a T he density is derived from the emission measure by assuming a source volume of ~ 8 X 10 26 cm 3 llKrucker et all 
l2010t) . 

b The kappa temperature k^T K is used for the kappa distribution fit. 

c The power-law index 8(= k + 0.5) is for the density distribution F(E) oc E~ s . 



Table 2. Fit result of the flatter component (colored blue in Figure [3]) 



fit parameters 


Model A 


Model B 


Model C 


Model C 




power-law 


thermal power-law 


kappa 


kappa 


emission measure EM, cm -3 




(1.6 ±0.3) X 10 46 


(8.6 ±6.2) x 10 46 


~ 8 x 10 46 


total density Ntot, cm -3 




(4.5 ±2.2) X 10 9 


(1.0 ±0.9) x 10 10 


~ 1 x 10 9 


temperature T, MK 




52 ±4 


28 ±9 


29 ± 18 


power-law index 8 


3.9 ±0.04 


3.8 ±0.2 


4.3 ±0.2 


4.2 ±0.4 



Note. — Same format as Table [T1 An error range (sigma level) is not shown when it exceeded the parameter value. 
The power-law index 8 is for the density distribution F(E) oc E~ s and can be converted from <5fd of the flux density 
distribution used in f_thin.pro (Models A & B; 8 = <5fd + 0.5) and from k of the number density distribution used 
in f_thin_kappa.pro (Models C & D; 8 = k + 0.5). 



- 19 - 




Fig. 4. — The ratio of flatter component (blue curves in Figure [3]) to steeper component 
(red curves in Figure [3]) derived by four different sets of models (histograms). The ratio is 
compared to the flux ratio (filled circles) derived from images (Figure [2]). We did not use 
the images of 14-17 and 17-20 keV ranges because of the unclear nature of the sources. 
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— Model A: power-law 

— Model B: thermal + power-law 

— Model C: kappa 



(a) 10 h i i | i i h 1 1 




0.4 

00:47:30 00:48:00 00:48:30 00:49:00 

31 December 2007 



Fig. 5. — Time variations of the fitted parameters for the flatter component (i.e. southern 
source) obtained by Model A (purple), Model B (green) and Model C (pink). From top to 
bottom are (a) the spectral index 5 and (b) the low-energy cutoff E c of the non-thermal tail 
of the southern source, (c) the emission measure EM and (d) temperature of the thermal 
core of the southern source, and (e) the reduced x 2 °f the model fits. The result by Model 
C is not shown because the error ranges (sigma levels) of EM and T exceeded the parameter 
values. While the peak spectra in Figure [3] were taken during the 8 s interval of 00:47:42- 
00:47:50 UT (indicated by the dashed lines), we fitted the data every 4s over the nearly 
two-minutes interval of 00:47:15 - 00:49:10 UT in this figure. 
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Table 3. Non-thermal fractions of electron number /energy densities in the southern source 



type of 


Model A 


Model B 


Model C 


density 


(power-law without thermal core) 


(power-law with thermal core) 


(kappa) 


number density 


44 % (>5 %) 


>2 % 


20 ± 1 % 


energy density 


88 % (>32 %) 


>12 % 


52+ 3 2 % 



Note. — See text for details of the estimations. For Model A, the lower limits of the number/energy 
densities (shown i n parenthesis) are b ased on the upper limit of the estimated number density of the 
ambient electrons llKrucker et alj feoiO). For Model B, the inequalities (>) are based on the fact that a 
range of values for the low-energy cutoff E c fit the data equally well and we used the cases of the highest 
E c . The result by Model C is not shown because the error ranges (sigma levels) of emission measure 
EM and temperature T exceeded the parameter values. 



